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CHARACTERIZATION OF SINGLE NUCLEOTIDE POLYMORPHISMS 

IN CODING REGIONS OF HUMAN GENES 

RELATED APPLICATION 

This application claims the benefit of U.S. Provisional Application Serial 
5 No. 60/127,248, filed March 31, 1999, the entire teachings of which are incorporated 
herein by reference. 

GOVERNMENT SUPPORT 

The invention was supported, in whole or in part, by grant 5-P50-HG00098- 
09 SNP from the National Institutes of Health (NCHGR) and grant 1-R01- 
10 HL61774-01 from the National Institutes of Health (NHLBI). The Government has 
certain rights in the invention. 



BACKGROUND OF THE INVENTION 

A major goal in human genetics is to understand the role of common genetic 
variants in susceptibility to common diseases (N. Risch and K. Merikangas, Science, 

15 275:1516 (1996.); E. S. Lander, Science, 274:536 (1996); F.S. Collins, et al. 
Science, 275:1580 (1997)). This requires assembling an extensive catalogue of 
single-nucleotide polymorphisms (SNPs) and performing systematic association 
studies for particular diseases. 

The human population has relatively limited genetic diversity, reflecting its 

20 young age and historically small size (F. J. Ayala et. al, Proc. Natl Acad. Set. 
97:6787 (1994)). Given the restricted nature of the allelic'spectrum, some authors 
have recently suggested that it should eventually be possible to collect all common 
SNPs in the human population and have hypothesized that such common variants 
may underlie much of the genetic risk of common disease (N. Risch and K. 

25 Merikangas, Science, 273:1516 (1996.); E. S. Lander, Science, 274:536 (1996); F.S. 
Collins, et al. Science, 278: 1 580 (1997)). This is in contrast to the situation for rare 
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Collins, et al, Science. 275:1580 (1997)). This is in contrast to the situation for rare 
genetic diseases, which are primarily caused by a large number of distinct alleles 
that are recent, rare and highly penetrant Important examples of associations to 
common (>1%) alleles include the ApoE4 allele in Alzheimer's disease, the Factor 
5 V u,dcn allele in deep-venous thrombosis, and the CCR5-A32 in resistance to HIV 
infection (A. M. Saunders et. at.. Neurology, 43:1467 (1993); R. M. Bertina, Nature. 
369:64 (1994); M. Dean et. al, Science, 275:1856 (1996)). The most relevant 
variants are likely to be those in coding and regulatory regions of genes. 



SUMMARY OF THE INVENTION 
1 0 As described herein, the nature of SNPs in the coding regions of human 

genes has been explored. SNPs were identified in 1 06 genes relevant to 
cardiovascular disease, endocrinology and neuropsychiatry, by screening an average 
of 1 14 independent alleles using two independent screening methods. To ensure 
high accuracy, all reported SNPs were confirmed by DNA sequencing. A total of 

15 545 SNPs were identified, including 395 coding-regions SNPs (cSNPs) divided 
roughly equally between those causing synonymous and non-synonymous changes. 
The cSNPs most likely to influence disease, those that alter the amino acid sequence 
of the encoded protein, show strikingly different properties: they occur at a lower 
rate and with lower allele frequencies. This likely reflects selection acting against 

20 deleterious alleles during human evolution. The lower allele frequency of cSNPs 
has important implications for the number of chromosomes that must be sampled to 
construct a comprehensive catalogue of human cSNPs. 

The invention relates to a gene which comprises a single nucleotide 
polymorphism at a specific location. In a particular embodiment the invention 

25 relates to the variant allele of a gene having a single nucleotide polymorphism, 
which variant allele differs from a reference allele by one nucleotide at the site(s) 
identified in Figures 5A-5QQQQQQQ. Complements of these nucleic acid 
segments are also included. The segments can be DNA or RNA, and can be double- 
or single-stranded. Segments can be, for example, 5-10, 5-15, 10-20, 5-25, 10-30, 

30 10-50 or 10-100 bases long. The invention further relates to gene products encoded 
by genes and oligonucleotides of the invention. 
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The invention further provides allele-specific oligonucleotides that hybridize 
to a gene comprising a single nucleotide polymorphism or to the complement of the 
gene. These oligonucleotides can be probes or primers. 

The invention further provides a method of analyzing a nucleic acid from an 
5 individual. The method determines which base is present at any one of the 

polymorphic sites shown in Figures 5A-5QQQQQQQ. Optionally, a set of bases 
occupying a set of the polymorphic sites shown in Figures 5A-5QQQQQQQ is 
determined. This type of analysis can be performed on a number of individuals, 
who are tested for the presence of a disease phenotype. The presence or absence of 
1 0 disease phenotype is then correlated with a base or set of bases present at the 
polymorphic site or sites in the individuals tested. 



BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a graph showing minor allele frequency by polymorphism type. 
The percentage of cSNPs having minor allele frequency classified as low (<5%), 
1 5 medium (5- 1 5%) or high (> 1 5%) frequency is displayed for synonymous, 
non-synonymous and non-coding SNPs. 

Figure 2 is a graph showing the distribution of nucleotide diversity. 
Normalized frequency of variant sites, ^, was calculated for the coding region of 
each gene. The graph shows the percentage of genes having & in the indicated 
20 range. 

Figures 3 A and 3B are a table showing a summary of polymorphisms in 106 
human genes described herein. Column 1 shows the name of the gene as used in 
Online Mendelian Inheritance in Man. Column 2 shows the number of coding base 
pairs screened. Column 3 shows the number of synonymous (or silent) 

25 polymorphisms identified. Column 4 shows the number of non-synonymous 

polymorphisms identified. Column 5 shows the number of non-coding base pairs 
screened. Column 6 shows the number of non-coding polymorphisms, including 
those in introns and untranslated regions (UTR), identified. 

Figure 4 is a table showing polymorphism rates for different classes of sites. 

30 Figures 5A-5QQQQQQQ are a table showing the specific polymorphisms 

identified in the genes studied as described herein. Column 1 shows the laboratory 
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designation for the polymorphism. Column 2 shows the name of the gene as used in 
Online Mendelian Inheritance in Man. Column 3 shows the reference nucleotide 
which occupies the polymorphic site in the reference allele. Column 4 shows the 
variant nucleotide which occupies the polymorphic site in the variant allele. 
5 Column 5 shows the reference amino acid encoded by the codon which contains the 
polymorphic site in the reference allele. Column 6 shows the variant or alternate 
amino acid encoded by the codon which contains the polymorphic site in the variant 
allele. Column 7 indicates whether the polymorphism is located in the coding or 
non-coding region of the gene. Column 8 shows the assay number in which the 
10 polymorphism was assessed. Columns 9 and 10 show the forward and reverse 
primers, respectively, which were used to identify the polymorphism. Column 1 1 
shows the sequence of the gene used in the assay, with the polymorphic site 
indicated by brackets and the primers shown in capital letters. Column 12 shows the 
total number of nucleotides given in Column 1 1 . 

1 5 DETAILED DESCRIPTION OF THE INVENTION 

There is a rich literature concerning nucleotide variation in model systems, 
particularly in Drosophila (E. N. Moriyama and J. R. Powell., MoL BioL EvoL, 
13:261 (1996)), but sequence variation in human genes has been studied only in 
limited ways. A small number of studies have focused on individual genes (such as 

20 beta-globin and lipoprotein lipase) in many individuals, and one study examined 49 
genes by comparing two independent sequences deposited in public databases (R. 
M. Harding et. al, Am. J. Hum. Genet, 60:112 (1997); D. A. Nickerson et. ai, 
Nature Genetics, 79:233 (1998); W. H. Li and L.A. Sadler, Genetics 129:513 
(1991)). To perform a more comprehensive survey, as described herein, a collection 

25 of 106 genes were selected whose protein products play important roles in the 

cardiovascular, endocrine and neurological systems (Figures 3A-3B and Figures 5A- 
5QQQQQQQ). Gene sequences were obtained from the Genbank and TIGR 
databases. Where multiple sequence depositions were available, a consensus 
sequence was derived. Determination of coding sequence, untranslated regions and 

30 intronic regions was based on annotation in the public database, although internal 
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checks were performed to ensure accurate determination of start and stop codons, 
open reading frames and the like. 

The genes were chosen because of their relevance to common, clinically 
significant diseases, such as coronary artery disease, diabetes, and schizophrenia. 
5 They encode proteins involved in coagulation, lipid metabolism, energy metabolism, 
neuroendocrine physiology, neurotransmission and central nervous system 
development. Variation in these genes was studied in a sample including 
Caucasians, African-Americans, African Pygmies and Asians, with an average of 
1 14 chromosomes screened for each gene. Of the samples screened, 30 were from 
10 Caucasian individuals, 14 from Asian, 10 African American and 7 Africans. The 
average number of individuals successfully screened for each gene was 57, with the 
precise number successfully screened varying among genes. Cell lines were 
obtained from Coriell Cell Repository, and DNA prepared according to standard 
protocols. In addition, 10 of the Caucasian samples used in this study were obtained 
1 5 as anonymous blood samples from the Physician's Health Study (gift of Charles 
Hennekens and J. Michael Gaziano). The sample size provides greater than 65% 
power to detect alleles with frequency of 1%. 

Overall, the sample of 1 14 chromosomes was screened for SNPs in a total of 
195.4 kb, consisting of 135.8 kb of coding regions and 59.6 kb from adjacent 
20 non-coding region (untranslated region (UTR) and introns). Sequences were 

amplified by the polymerase chain reaction (PCR) and screened by two independent 
methods. The first method involved hybridization of labeled PCR products to 
variant detector arrays (VDAs) (that is, high density DNA probe arrays containing 
oligonucleotides specific for the sequences under study) (M. Chee et. al. Science. 
25 274:610 (1996); D. G. Wang et, al., Science, 280:1071 (1998)); variant sequences 
typically give rise to altered hybridization patterns. These chips contained variant 
detector arrays (VDA) (M. Chee et aU Science 274:610 (1996)). 

Using VDAs, candidate SNPs were identified using a combination of three 
algorithms followed by visual inspection. For each base position and strand queried 
30 there are four VDA features: one contains the expected base (the reference 
sequence) in the central position and the other three features contain central 
substitution bases (in the background of the reference sequence). The base-calling 
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algorithm looked for positions at which hybridization to a substitution base gives a 
stronger signal than the reference base. The second algorithm (mutant fraction) 
examined the reference base and each one of the substitution bases in turn and 
calculates the fraction of signal present in the non-reference base. The final 
5 algorithm (footprint detection) depends upon a loss of signal at the reference 

positions surrounding a nucleotide substitution. These algorithms are combined to 
yield a confidence score of "certain" or "likely" for each candidate polymorphism. 
Two analysts independently scored the data, and candidate polymorphisms found by 
either observer were included in subsequent confirmation tests. PCR assays 
1 0 spanning each exon were designed using Primer 3.0 release 0.7. PCR was 

performed according to standard protocols, and assays destined to be hybridized to 
the same chip design were pooled together. Chip samples were prepared and 
hybridized as described in D.G. Wang et al. {Science 250:1077 (1998)), except that 
pools consisting of about 100 assays contained 5-6 Jig of amplified material. In all, 
15 854 assays (average size of 300 bp, covering 106 genes) were amplified from each 
individual and were hybridized to 12 distinct chip designs. The probe arTays were 
designed to query only the coding sequence for some genes, while other genes 
contained the entire mRNA and/or surrounding intron (Figures 3A-3B). The second 
method involved subjecting PCR products to Denaturing HPLC (dHPLC) (P. J. 
20 Oefiier and P. A. Underhill, Am. J. Hum. Genet, J7.A266 (1995)) at a critical 

temperature; heterozygous individuals typically give rise to heteroduplex products 
with altered denaturation and migration properties. 

Sequences were amplified as above except that the final extension in the 
PCR protocol was followed by denaturation and slow reannealing to allow 
25 heteroduplex formation. A total of 6 |JLl of each individual PCR product was 

injected into Wave DNA Fragment Analysis System (Transgenomic). A total of 592 
of the VDA assays (covering the 89 genes attempted with this method) were 
successfully screened by DHPLC. Only assays of >160 base pairs were used for 
DHPLC, because shorter assays performed unreliably for mutation detection. The 
30 DHPLC parameters (percentage of acetonitrile, column temperature) used for each 
fragment were automatically calculated using a novel predictive algorithm, and 
DHPLC traces were analyzed using the clustering program ASH v2.0. A scoring 
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algorithm was developed based upon the similarity score by ASHv2.0 and contour 
of the elution profile. 

Because both screening methods can generate to a significant number of 
false positives, it was important to confirm every reported SNP. Samples implicated 
5 by either method as containing a candidate SNP were thus subjected to fluorescent 
dideoxy sequencing, either to confirm the presence of the SNP (in the case of the 
chip) or to identify and confirm the presence of the SNP (in the case of DHPLC). 
Such confirmation proved essential for eliminating false positives. 

Candidate SNPs were either validated (if found by VDAs) or identified (if 

10 implicated by DHPLC) by DNA sequencing. For this purpose, sequences were 

amplified with PCR primers tailed with standard Ml 3 sequencing sites (-21 forward 
and -28 reverse) and conventional dye-primer sequencing was performed on ABI 
377 sequencers. For candidate SNPs discovered by VDAs, one individual was 
chosen (a candidate homozygous variant, when available, or a candidate 

15 heterozygote) and sequencing was performed on one strand to confirm by visual 
inspection the presence of the SNP at the indicated position. For ampiicons found 
to be polymorphic by DHLPC, two individuals were selected representing each 
distinct elution pattern observed and were sequenced on both strands to discover the 
variant base or bases. Sequences were base-called by the Phred program, assembled 

20 by the Phrap program, and polymorphism candidates were identified by the 

PolyPhred program (D. A. Nickerson eu aL, NAR, 25:2145 (1997). All results were 
visually inspected by at least two observers. 

The overall false positive rate for VDAs was 45%. The rate was much lower 
(about 10%) for certain chip designs, synthesis protocols, and for candidate 

25 polymorphisms scored as "certain." The false positive rate among fragments 
displaying an altered elution pattern by DHPLC was similar (40%). The false 
positive rates reflect the thresholds employed for declaring a candidate SNP, which 
were chosen to ensure high sensitivity. 

A total of 545 SNPs were identified in the 195 kb surveyed, consisting of 

30 150 non-coding SNPs and 395 cSNPs, Results from these studies are shown in the 
Figures. The complete data are available on the web site 

http://www.genome.wi.mit.edu/cvar_snps; access to this website can be gained 
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using the guestname "snp_pilot" and the password "noynek". In the future, access to 
this website may be available to the public, and thus, no guestname or password may 
be needed. 

To directly determine the false-negative rate of the screen, conventional 
5 DNA sequencing was performed on ten of the genes (THPO, TBAX2R, PTHLH, 
IGF2, HTR2A, HTR1A, GHR, GABRB1, F10, and CYP1 1B1) spanning 25.2 kb in 
twenty individuals. Sequencing was performed on both strands using dye-primer 
chemistry and sequence traces were interpreted using PolyPhred (D.A. Nickerson et. 
aL, NAR, 25:2145 (1997)). VDA analysis identified 85% of variants found by direct 
10 sequencing , while DHPLC identified 87% of the variants found by direct 

sequencing. In regions screened by both VDAs and DHPLC, the combination of the 
two methods identified 100% of the polymorphisms found by direct sequencing. 

Overall, about one-third of individuals were screened with both methods, and 
one-third were screened with each of the two methods alone. (For some genes, the 
15 non-coding regions were screened only by DHPLC.) It is estimated that the false 
negative rate over the entire study to be about 15% for regions screened by one 
method, and negligible for sequences screened by both methods. The total number 
of true polymorphisms not identified is estimated to be less than 10%. 

A SNP survey can be characterized in terms of either K, the observed 
20 number of variant sites, or p, the observed heterozygosity per bp. Because K 

increases with the number of chromosomes (n) studied and the total sequence length 
L, it is preferable to use the normalized number of variant sites 

© = K/ i" 1 ) L which corrects for sample size. Under the neutral 

theory of molecular evolution and infinite sites model, 0 and n are both estimators 
25 of the population genetic parameter 0 = 4N\l (Li, Molecular Evolution, Sinauer 

Associates (1997), Canada). 

SNPs were found at a similar overall frequency in coding and non-coding 

regions. SNPs in coding region occurred at a frequency of 1 per 344 bp, 

corresponding to & = 5.47 x 10" 4 and tz = 5.07 x 10" 4 . Interestingly, SNPs were 
30 observed in non-coding DNA at a similar frequency of 1 per 397 bp. The 



WO 00/58519 



-9- 



PCT/US00/08440 



normalized number of variant sites was 6 = 4.93 x 10\ and the mean 
heterozygosity (it) = 5.05 x IV 4 (Figure 4). Calculations of * involve allele 
frequencies. Polymorphisms identified by DHPLC alone were excluded because we 
did not sequence all of the samples showing a variant DHPLC pattern and thus could 

5 not be certain of allele frequency. The estimates of TZ were thus based on 411 of 545 
polymorphisms. Although the VDAs were designed for polymorphism discovery 
rather than genotyping, the estimated allele frequencies proved to be quite accurate. 
Specifically, genotyping assays (employing single-base extension assays) for 25 
SNPs yielded allele frequencies that differed by an average of only 2% from those 

10 estimated on the basis of genotypes inferred from the VDA. For both classes, the 
similar values for 6 and it is consistent with a population evolving according to 

neutral expectations. 

The 395 cSNPs were roughly equally divided between synonymous (203 
cSNPs) and non-synonymous (192 cSNPs) changes. Since approximately two-thirds 
15 of random mutations would alter an amino acid, the fact that non-synonymous 

cSNPs comprise slightly less than half of the cSNPs implies strong selection against 
amino-acid altering changes. To address this issue more directly, the nucleotide 
diversity was examined at four-fold degenerate sites, two-fold degenerate sites, and 
non-degenerate sites. Changes at four-fold degenerate sites produce only 
20 synonymous changes, while those at non-degenerate sites are always 

non-synonymous. Nucleotide diversity (*> was 9.64 x 10"< at four-fold degenerate 
sites, 6.85 x Iff 4 at two-fold degenerate sites, and 3.70 x 10" 4 at non-degenerate sues. 
Assuming that mutations occur at an equal rate at both classes of sites, 
non-synonymous variants survive to be detected in such a survey at only 38% of the 

25 rate of synonymous changes. 

The force of selection is also evident in comparing non-synonymous cSNPs 
causing a non-conservative amino acid alteration with those causing a conservative 
amino-acid change. Conservative and non-conservative amino acid substitutions 
were defined for this analysis according to the BLOSUM62 matrix, used in sequence 

30 comparison(S.HenikoffandJ.G.Henikoff,P^5.S9:10915(1992). Conservative 

changes were those having a positive or neutral sign in the matrix, while 
non-conservative changes were those having a negative value Non-conservative 
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cSNPs represent only 36% of the non-synonymous cSNPs, whereas randomly 
distributed mutations would be expected to produce a higher proportion (52%) of 
non-conservative changes. The proportion of non-synonymous SNPs expected to 
cause a non-conservative amino acid substitution was determined based on the 
5 actual codon usage in the 1 06 genes studied, the known frequencies of transitions 
and transversions, and the definition of non-conservative changes employed in the 
BLOSUM62 matrix. This implies that non-conservative cSNPs survive to be 
detected in such a survey at only about half of the rate of conservative, 

non-synonymous cSNPs. 

1 o The various types of SNPs differ not only in the rate of their occurrence, but 

also in the frequency of their minor alleles. This can be seen in several ways. When 
SNPs are classified according to whether the frequency of the minor allele was high 
(>15%),intermediate (5-15%) or low (<5%), it is clear that the non-synonymous 
cSNPs were enriched in low frequency alleles compared to the rest of the collection 

15 (Figure 1). The distribution of non-synonymous allele frequencies was significantly 
different than that of synonymous changes (p=0.02, Kolmogorov-Smirnov test). 
Indeed, more than half (58%) of non-synonymous cSNPs were found at a frequency 
below 5%, with this effect evident for both conservative and non-conservative 
substitutions. 

20 The effect of selection can also be inferred by considering the average 

frequency of the minor allele: it is 8% for non-conservative cSNPs, 1 1% for 
conservative but non-synonymous cSNPs, and 14% for both synonymous cSNPs 
and non-coding SNPs. In addition, the lower allele frequency of non-synonymous 
cSNPs is reflected in the fact that the heterozygosity Ti is lower than the normalized 

25 rate of variant sites £ for this class of SNPs (Figure 4). This divergence is in the 
direction predicted by the action of purifying selection, although it falls short of 
statistical significance. Tajima's D was non-significant. (F. Tajima, Genetics. 
723:545 (1989). 

The distribution of SNPs among the 106 genes was explored, with an eye 
30 toward detecting differential effects of selection among genes. The number of 
cSNPs per gene ranged from 37 for Factor V to 0 for thirteen of the genes, and the 
normalized rate, B, similarly showed considerable variation (Figure 2). The 
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observed variation in nucleotide diversity is similar in magnitude to that observed 
for Drosophila (E. N. Moriyama and J. R. Powell., Mol Biol Evol, 75:261 (1996)). 
Variation among genes could be due to many factors (D. J. Begun and C. F. 
Aquadro, Nature, 356:519 (1993); Nachmaner. ai, Genetics, 750:1133 (1998)). 
5 The fact that non-synonymous cSNPs show a somewhat wider variation than 

synonymous cSNPs (the coefficient of variation is 20% larger for the former class) 
is consistent with differences in selective constraints among loci, but the difference 
falls well below statistical significance. A variety of population genetic tests are 
available for testing selection at individual loci (M. L. Wayne and K. L. Simonson, 

10 Trends and Ecology and Evolution, 73:236 (1998)). 

The age of a SNP allele has important implications for its use in human 
genetic studies. Recently-occurring SNP alleles are more likely to show extensive 
linkage disequilibrium (retention of the ancestral haplotype on which they arose) as 
compared to older SNPs. Such linkage disequilibrium can provide a powerful tool 

15 in identifying disease genes (E. S. Lander, N.J. Schork, Science, 265:2037 (1994)). 
Although the precise age of the SNPs could not be assessed from these studies, 
characterization of which allele preceded human speciation and which arose 
thereafter was sought. To determine the ancestral human allele, each corresponding 
gene was sequenced from the common chimpanzee (P. troglodytes). Each assay 

20 used in the human survey was amplified from a single chimpanzee (DNA gift of 
Kristin Ardlie) and subjected to dye-primer sequencing on both strands. A single 
chimpanzee sample will accurately reveal the ancestral allele except in cases where 
the site has mutated and fixed during the chimpanzee evolution or is polymorphic in 
the chimpanzee population and happened to be homozygous for the non-ancestral 

25 allele. These two cases are quite rare (probably less than 2%) and thus have been 

neglected for the purpose of estimating overall rates. A human allele was considered 
to be ancestral if it was present in the homozygous state in the chimpanzee sample. 
A total of 136 kb of chimpanzee sequence was obtained, revealing an inter-species 
divergence of 0.6% in the regions studied. 

30 An elegant result in theoretical population genetics predicts that the 

probability that a neutral allele represents the ancestral state should be equal to its 
frequency in the population (G. A. Watterson and H. A. Guess, Theoretical 
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Population Biology, 77:141 (1977)). The minor allele should thus represent the 
ancestral state in a predictable proportion of cases. The ancestral allele and minor 
allele frequency was determined for 267 of the reported SNPs. For 3 of the 267 
SNPs, the chimpanzee was homozygous for a third allele differing from both of the 
5 current human alleles. This is consistent with the overall 0.6% nucleotide sequence 
divergences seen between human and chimpanzee. Among polymorphisms with a 
minor allele frequency below 10%, the average allele frequency was 3% and the 
proportion that was ancestral was 7% (1 1/158) of cases. Among polymorphisms 
with minor alleles exceeding 10%, the mean frequency was 28% and the proportion 

1 0 that were ancestral was 32% (35/1 09). These results thus agree remarkably well 
with the theoretical prediction, providing the first reported test of this prediction in 
humans. It therefore follows that the minor SNP allele need not be the younger 
allele; this has implications for linkage disequilibrium mapping. 

The distribution of SNPs among Caucasian, African- American, African and 

15 Asian samples was also examined. Although the vast majority of SNPs were seen in 
multiple groups, there was a statistically significant excess of SNPs that were seen in 
only one of the sub-groups. The probability that a SNP occurring k>l times in an 
overall sample of n individuals would be found entirely within a given subset of m 
individuals is B(n,k)/B(m,k), where B(x,y) is the binomial coefficient x!/(x-y)!y!. In 

20 this fashion, the probability that each individual SNP would be confined to a 

particular ethnic subgroup within the sample was calculated and these probabilities 
were summed to obtain the number of SNPs expected to be confined to the group 
within the sample. The fact that a SNP is found only within one group in the sample 
does not necessarily imply that it is private to that group within the general 

25 population, owing to the small sample size, but it can be used as an indication of 
substructure. The number of SNPs with k > 1 confined to the, African- Americans, 
African Pygmies, Caucasians, and Asians was 17, 17, 12, and 9, as compared to 
expectations of 3.02, 1.34 , 8.62, and 1.81. Not surprisingly, the greatest excess was 
seen for SNPs found in the African-American and African samples. The presence of 

30 population substructure implies that construction of a comprehensive SNP database 
should employ a diverse set of DNA samples. 
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The results of this survey provide a fundamental description of sequence 
variation in the coding regions of human genes. These data indicate that two copies 
of a gene chosen from the human population will differ by roughly one base in 2 kb, 
corresponding to somewhat less than one heterozygous base within the coding 
5 region of a typical gene. In general, there are only a handful of such cSNPs per gene 
that exhibit allele frequencies of at least a few percent. Accounting for both the 
different rate and frequency of non-synonymous SNPs, only about 40% of these 
observed changes will alter the encoded amino acid. The action of purifying 
selection during human evolution is evident from the comparatively lower rate of 
10 non-synonymous cSNPs, and especially of those that create a non-conservative 
change. It is clear that non-synonymous cSNPs not only occur less often, but also 
have lower minor allele frequencies: 60% of non-synonymous cSNPs, the class 
likely to have the most dramatic effects on proteins, display a minor allele frequency 
below 5%. 

1 5 The relative rarity of cSNPs has important implications for efforts to produce 

large catalogues of human variants. It has been proposed that most human SNPs 
could be found by performing shotgun sequencing on a handful of individuals (J. L. 
Weber and E. W. Myers, Genome Research, 7:401 (1997); J. C. Venter et. aL, 
Science, 280:1540 (1998)). Although such a project will surely identify many SNPs, 

20 results described herein suggest that the small sample size will likely fail to identify 
the vast majority of cSNPs likely to have the most important biological 
consequences, owing to their lower average allele frequencies. A comprehensive 
collection of the common, non-conservative cSNPs may require surveying 50-100 
chromosomes. Because coding sequence represents only about 3% of the genome, it 

25 may prove inefficient to obtain such deep coverage of cSNPs by shotgun sequencing 
of genomic DNA. Instead, it may be more efficient to perform shotgun sequencing 
on cDNA libraries from multiple individuals or to amplify genes from multiple 
individuals, as done here. 

Interestingly, a similar rate of polymorphism in coding and non-coding DNA 

30 was found. Furthermore, the observed rate of nucleotide diversity at four- fold 
degenerate sites was nearly twice that in adjacent non-coding regions, and over 
twice that at non-degenerate sites (Figure 4). Similar results have been reported for 
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Drosophila (E. N. Moriyama and J. R. Powell., MoL BioL EvoL, 13:261 (1996)) and 
for a smaller human data set by Li and Sadler (R. M. Harding et. al„ Am. J, Hum. 
Genet., 60:112 (1997); D. A. Nickerson et. al , Nature Genetics, 19:233 (1998); 
W. H. Li and L.A. Sadler, Genetics 129:513 (1991)), who observed over three times 
5 the nucleotide diversity at four-fold degenerate sites (& = 1 1 x 10" 4 ), as compared to 
that in both untranslated regions and non-degenerate sites (^ = 3x 10" 4 ). These 
observations suggest that non-coding DNA adjacent to coding regions may be 
functionally constrained to a surprising degree. 

SNPs can be used to search for genes underlying complex traits in two 
10 distinct ways: linkage disequilibrium (LD) studies and association studies (E. S. 
Lander, N.J. Schork, Science, 265:2031 (1994)). Genome-wide LD studies involve 
using a dense collection of SNPs as markers to search for an ancestral haplotype 
carrying a disease-susceptibility allele. Such studies cannot be undertaken without 
the availability of an extremely dense SNP map and their potential for success 
1 5 depends sensitively on many population genetic assumptions. Association studies 
are more straightforward because they directly test the hypothesis that a specific 
SNP increases disease risk. They make few assumptions, and require only the 
availability of a suitable database of appropriate SNPs. In the near term, focusing on 
cSNPs is likely to be most productive inasmuch as the class is easily recognized (in 
20 contrast to regulatory polymorphisms) and is likely to contain a significant 
proportion of the disease-susceptibility alleles. 

The present invention relates to a gene which comprises a single nucleotide 
polymorphism (SNP) at a specific location. The gene which includes the SNP has at 
least two alleles, referred to herein as the reference allele and the variant allele. The 
25 reference allele (prototypical or wild type allele) has been designated arbitrarily and 
typically corresponds to the nucleotide sequence of the gene which has been 
deposited with GenBank or TIGR under a given Accession number. The variant 
allele differs from the reference allele by one at least one nucleotide at the site(s) 
identified in Figures 5A-5QQQQQQQ. The present invention also relates to variant 
30 alleles of the described genes and to complements of the variant alleles. The 
invention further relates to portions of the variant alleles and portions of 
complements of the variant alleles which comprise (encompass) the site of the SNP 



WO 00/58519 



-15- 



PCT/US00/08440 



and are at least 5 nucleotides in length. Portions can be, for example, 5-10, 5-15, 10- 
20, 5-25, 10-30, 10-50 or 10-100 bases long. For example, a portion of a variant 
allele which is 5 nucleotides in length includes the single nucleotide polymorphism 
(the nucleotide which differs from the reference allele at that site) and four 
5 additional nucleotides which flank the site in the variant allele. These nucleotides 
can be on one or both sides of the polymorphism. Polymorphisms which are the 
subject of this invention are defined in Figures 5A-5QQQQQQQ with respect to the 
reference sequence deposited in GenBank under the Accession number indicated. 
For example, the invention relates to a portion of a gene (e.g., AADC) having a 

10 partial nucleotide sequence as shown in Figures 5A-5QQQQQQQ comprising a 
single nucleotide polymorphism at a specific position. The reference nucleotide for 
AADC is shown in column 3 and the variant nucleotide is shown in column 4 of 
Figures 5A-5QQQQQQQ. The nucleotide sequences of the invention can be 
double- or single-stranded. 

15 The invention further provides allele-specific oligonucleotides that hybridize 

to a gene comprising a single nucleotide polymorphism or to the complement of the 
gene. These oligonucleotides can be probes or primers. 

The invention further provides a method of analyzing a nucleic acid from an 
individual. The method determines which base is present at any one of the 

20 polymorphic sites shown in Figures 5A-5QQQQQQQ. Optionally, a set of bases 
occupying a set of the polymorphic sites shown in Figures 5A-5QQQQQQQ is 
determined. This type of analysis can be performed on a number of individuals, 
who are tested for the presence of a disease phenotype. The presence or absence of 
disease phenotype is then correlated with a base or set of bases present at the 

25 polymorphic site or sites in the individuals tested. 

An oligonucleotide of this invention can be DNA or RNA, and single- or 
double-stranded. Oligonucleotides can be naturally occurring or synthetic, but are 
typically prepared by synthetic means. Preferred oligonucleotides of the invention 
include segments of DNA, or their complements, which include any one of the 

30 polymorphic sites shown in Figures 5A-5QQQQQQQ. The segments can be 

between 5 and 250 bases, and, in specific embodiments, are between 5-10, 5-20, 10- 
20, 10-50, 20-50 or 10-100 bases. The polymorphic site can occur within any 
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position of the segment. The segments can be from any of the allelic forms of DNA 
shown in Figures 5A-5QQQQQQQ. 

As used herein, the terms "nucleotide" and "nucleic acid" are intended to be 
equivalent. The terms "nucleotide sequence", "nucleic acid sequence" "nucleic acid 
5 molecule" and "segment" are intended to be equivalent. 

Hybridization probes are oligonucleotides which bind in a base-specific 
manner to a complementary strand of nucleic acid. Such probes include peptide 
nucleic acids, as described in Nielsen et aL, Science 254, 1497-1500 (1991). Probes 
can be any length suitable for specific hybridization to the target nucleic acid 

10 sequence. The most appropriate length of the probe may vary depending upon the 
hybridization method in which it is being used; for example, particular lengths may 
be more appropriate for use in microfabricated arrays, while other lengths may be 
more suitable for use in classical hybridization methods. Suitable probes and 
primers can range from about 5 nucleotides to about 30 nucleotides in length. For 

15 example, probes and primers can be 5, 6, 8, 10, 12, 14, 16, 18, 20, 22, 24, 25, 26, 28 
or 30 nucleotides in length. The probe or primer preferably contains at least one 
polymorphic site occupied by any of the possible variant nucleotides. The 
nucleotide sequence can correspond to the coding sequence of the allele or to the 
complement of the coding sequence of the allele. 

20 As used herein, the term "primer" refers to a single-stranded oligonucleotide 

which acts as a point of initiation of template-directed DNA synthesis under 
appropriate conditions (e.g., in the presence of four different nucleoside 
triphosphates and an agent for polymerization, such as, DNA or RNA polymerase or 
reverse transcriptase) in an appropriate buffer and at a suitable temperature. The 

25 appropriate length of a primer depends on the intended use of the primer, but 
typically ranges from 15 to 30 nucleotides. Short primer molecules generally 
require cooler temperatures to form sufficiently stable hybrid complexes with the 
template. A primer need not reflect the exact sequence of the template, but must be 
sufficiently complementary to hybridize with a template. The term primer site refers 

30 to the area of the target DNA to which a primer hybridizes. The term primer pair 
refers to a set of primers including a 5' (upstream) primer that hybridizes with the 5' 



WO 00/58519 



-17- 



PCT/USOO/08440 



end of the DNA sequence to be amplified and a 3' (downstream) primer that 
hybridizes with the complement of the 3* end of the sequence to be amplified. 

As used herein, linkage describes the tendency of genes, alleles, loci or 
genetic markers to be inherited together as a result of their location on the same 
5 chromosome. It can be measured by percent recombination between the two genes, 
alleles, loci or genetic markers. 

As used herein, polymorphism refers to the occurrence of two or more 
genetically determined alternative sequences or alleles in a population. A 
polymorphic marker or site is the locus at which divergence occurs. Preferred 

10 markers have at least two alleles, each occurring at frequency of greater than 1%, 
and more preferably greater than 10% or 20% of a selected population. A 
polymorphic locus may be as small as one base pair. Polymorphic markers include 
restriction fragment length polymorphisms, variable number of tandem repeats 
(VNTR's), hypervariable regions, minisatellites, dinucleotide repeats, trinucleotide 

15 repeats, tetranucleotide repeats, simple sequence repeats, and insertion elements 

such as Alu. The first identified allelic form is arbitrarily designated as the reference 
form and other allelic forms are designated as alternative or variant alleles. The 
allelic form occurring most frequently in a selected population is sometimes referred 
to as the wildtype form. Diploid organisms may be homozygous or heterozygous 

20 for allelic forms. A diallelic or biallelic polymorphism has two forms. A triallelic 
polymorphism has three forms. 

By altering amino acid sequence, SNPs may alter the function of the encoded 
proteins. The discovery of the SNP facilitates biochemical analysis of the variants 
and the development of assays to characterize the variants and to screen for 

25 pharmaceutical that would interact directly with on or another form of the protein. 
SNPs (including silent SNPs) may also alter the regulation of the gene at the 
transcriptional or post-transcriptional level. SNPs (including silent SNPs) also 
enable the development of specific DNA, RNA, or protein-based diagnostics that 
detect the presence or absence of the polymorphism in particular conditions. 

30 A single nucleotide polymorphism occurs at a polymorphic site occupied by 

a single nucleotide, which is the site of variation between allelic sequences. The site 
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is usually preceded by and followed by highly conserved sequences of the allele 
(e.g., sequences that vary in less than 1/100 or 1/1000 members of the populations). 

A single nucleotide polymorphism usually arises due to substitution of one 
nucleotide for another at the polymorphic site. A transition is the replacement of 
5 one purine by another purine or one pyriraidine by another pyrimidine. A 

transversion is the replacement of a purine by a pyrimidine or vice versa. Single 
nucleotide polymorphisms can also arise from a deletion of a nucleotide or an 
insertion of a nucleotide relative to a reference allele. Typically the polymorphic 
site is occupied by a base other than the reference base. For example, where the 

10 reference allele contains the base "T M at the polymorphic site, the altered allele can 
contain a "C\ "G" or "A" at the polymorphic site. 

Hybridizations are usually performed under stringent conditions, for 
example, at a salt concentration of no more than 1 M and a temperature of at least 
25°C. For example, conditions of 5X SSPE (750 mM NaCl, 50 mM NaPhosphate, 

15 5 mM EDTA, pH 7.4) and a temperature of 25-30°C, or equivalent conditions, are 
suitable for allele-specific probe hybridizations. Equivalent conditions can be 
determined by varying one or more of the parameters given as an example, as known 
in the art, while maintaining a similar degree of identity or similarity between the 
target nucleotide sequence and the primer or probe used. 

20 The term "isolated" is used herein to indicate that the material in question 

exists in a physical milieu distinct from that in which it occurs in nature. For 
example, an isolated nucleic acid of the invention may be substantially isolated with 
respect to the complex cellular milieu in which it naturally occurs. In some 
instances, the isolated material will form part of a composition (for example, a crude 

25 extract containing other substances), buffer system or reagent mix. In other 

circumstance, the material may be purified to essential homogeneity, for example as 
determined by PAGE or column chromatography such as HPLC. Preferably, an 
isolated nucleic acid comprises at least about 50, 80 or 90 percent (on a molar basis) 
of all macromolecular species present. 
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I. Analysis of Polymorphisms 

A. Preparation of Samples 

Polymorphisms are detected in a target nucleic acid from an individual being 
analyzed. For assay of genomic DNA, virtually any biological sample (other than 
5 pure red blood cells) is suitable. For example, convenient tissue samples include 
whole blood, semen, saliva, tears, urine, fecal material, sweat, buccal, skin and hair. 
For assay of cDNA or mRNA, the tissue sample must be obtained from an organ in 
which the target nucleic acid is expressed. For example, if the target nucleic acid is 
a cytochrome P450, the liver is a suitable source. 

10 Many of the methods described below require amplification of DNA from 

target samples. This can be accomplished by e.g., PCR. See generally PCR 
Technology: Principles and Applications for DNA Amplification (ed. H.A. Erlich, 
Freeman Press, NY, NY, 1992); PCR Protocols: A Guide to Methods and 
Applications (eds. Innis, et al., Academic Press, San Diego, CA, 1990); Mattila et 

1 5 aL, Nucleic Acids Res. 1 9, 4967 ( 1 99 1 ); Eckert et aL % PCR Methods and 

Applications 1, 17 (1991); PCR (eds. McPherson et aL f IRL Press, Oxford); and 
U.S. Patent 4,683,202. 

Other suitable amplification methods include the ligase chain reaction (LCR) 
(see Wu and Wallace, Genomics 4, 560 (1989), Landegren et aL, Science 241, 1077 

20 (1988), transcription amplification (Kwoh et al., Proc. Natl. Acad. Sci. USA 86, 
1 173 (1989)), and self-sustained sequence replication (Guatelli et al., Proc. Nat. 
Acad. Sci. USA, 87, 1874 (1990)) and nucleic acid based sequence amplification 
(NASB A). The latter two amplification methods involve isothermal reactions based 
on isothermal transcription, which produce both single stranded RNA (ssRNA) and 

25 double stranded DNA (dsDNA) as the amplification products in a ratio of about 30 
or 100 to 1, respectively. 

B. Detection of Polymorphisms in Target DNA 

There are two distinct types of analysis of target DNA for detecting 
polymorphisms. The first type of analysis, sometimes referred to as de novo 
30 characterization, is carried out to identify polymorphic sites not previously 

characterized (i.e., to identify new polymorphisms). This analysis compares target 
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sequences in different individuals to identify points of variation, i.e., polymorphic 
sites. By analyzing groups of individuals representing the greatest ethnic diversity 
among humans and greatest breed and species variety in plants and animals, patterns 
characteristic of the most common alleles/haplotypes of the locus can be identified, 
5 and the frequencies of such alleles/haplotypes in the population can be determined. 
Additional allelic frequencies can be determined for subpopulations characterized by 
criteria such as geography, race, or gender. The de novo identification of 
polymorphisms of the invention is described in the Examples section. The second 
type of analysis determines which form(s) of a characterized (known) polymorphism 
10 are present in individuals under test. There are a variety of suitable procedures, 
which are discussed in turn. 

1 . Allele-Specific Probes 

The design and use of allele-specific probes for analyzing polymorphisms is 
described by e.g., Saiki et aL, Nature 324, 163-166 (1986); Dattagupta, EP 235,726, 

15 Saiki, WO 89/1 1548. Allele-specific probes can be designed that hybridize to a 
segment of target DNA from one individual but do not hybridize to the 
corresponding segment from another individual due to the presence of different 
polymorphic forms in the respective segments from the two individuals. 
Hybridization conditions should be sufficiently stringent that there is a significant 

20 difference in hybridization intensity between alleles, and preferably an essentially 
binary response, whereby a probe hybridizes to only one of the alleles. Some probes 
are designed to hybridize to a segment of target DNA such that the polymorphic site 
aligns with a central position (e.g., in a 15-mer at the 7 position; in a 16-mer, at 
either the 8 or 9 position) of the probe. This design of probe achieves good 

25 discrimination in hybridization between different allelic forms. 

Allele-specific probes are often used in pairs, one member of a pair showing 
a perfect match to a reference form of a target sequence and the other member 
showing a perfect match to a variant form. Several pairs of probes can then be 
immobilized on the same support for simultaneous analysis of multiple 

30 polymorphisms within the same target sequence. 
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2. Tiling Arrays 

The polymorphisms can also be identified by hybridization to nucleic acid 
arrays, some examples of which are described in WO 95/1 1995. One form of such 
arrays is described in the Examples section in connection with de novo identification 
5 of polymorphisms. The same array or a different array can be used for analysis of 
characterized polymorphisms. WO 95/1 1995 also describes subarrays that are 
optimized for detection of a variant form of a precharacterized polymorphism. Such 
a subarray contains probes designed to be complementary to a second reference 
sequence, which is an allelic variant of the first reference sequence. The second 

10 group of probes is designed by the same principles as described in the Examples, 
except that the probes exhibit complementarity to the second reference sequence. 
The inclusion of a second group (or further groups) can be particularly useful for 
analyzing short subsequences of the primary reference sequence in which multiple 
mutations are expected to occur within a short distance commensurate with the 

15 length of the probes (e.g., two or more mutations within 9 to 21 bases). 

3. Allele-Specific Primers 

An allele-specific primer hybridizes to a site on target DNA overlapping a 
polymorphism and only primes amplification of an allelic form to which the primer 
exhibits perfect complementarity. See Gibbs, Nucleic Acid Res. 1 7, 2427-2448 

20 (1989). This primer is used in conjunction with a second primer which hybridizes at 
a distal site. Amplification proceeds from the two primers, resulting in a detectable 
product which indicates the particular allelic form is present. A control is usually 
performed with a second pair of primers, one of which shows a single base 
mismatch at the polymorphic site and the other of which exhibits perfect 

25 complementarity to a distal site. The single-base mismatch prevents amplification 
and no detectable product is formed. The method works best when the mismatch is 
included in the 3'-most position of the oligonucleotide aligned with the 
polymorphism because this position is most destabilizing to elongation from the 
primer (see, e.g., WO 93/22456). 
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4. Direct-Sequencing 

The direct analysis of the sequence of polymorphisms of the present 
invention can be accomplished using either the dideoxy chain termination method or 
the Maxam Gilbert method (see Sambrook et al. 9 Molecular Cloning, A Laboratory 
5 Manual (2nd Ed., CSHP, New York 1 989); Zyskind et al, Recombinant DNA 
Laboratory Manual (Acad. Press, 1988)). 

5. Denaturing Gradient Gel Electrophoresis 

Amplification products generated using the polymerase chain reaction can be 
analyzed by the use of denaturing gradient gel electrophoresis. Different alleles can 
1 0 be identified based on the different sequence-dependent melting properties and 
electrophoretic migration of DNA in solution. Erlich, ed., PCR Technology, 
Principles and Applications for DNA Amplification, (W.H. Freeman and Co, New 
York, 1992), Chapter 7. 

6. Single-Strand Conformation Polymorphism Analysis 

15 Alleles of target sequences can be differentiated using single-strand 

conformation polymorphism analysis, which identifies base differences by alteration 
in electrophoretic migration of single stranded PCR products, as described in Orita 
et al. f Proc. Nat. Acad. Sci. 86, 2766-2770 (1989). Amplified PCR products can be 
generated as described above, and heated or otherwise denatured, to form single 

20 stranded amplification products. Single-stranded nucleic acids may refold or form 
secondary structures which are partially dependent on the base sequence. The 
different electrophoretic mobilities of single-stranded amplification products can be 
related to base-sequence differences between alleles of target sequences. 

II. Methods of Use 

25 After determining polymorphic form(s) present in an individual at one or 

more polymorphic sites, this information can be used in a number of methods. 
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A. Forensics 

Determination of which polymorphic forms occupy a set of polymorphic 
sites in an individual identifies a set of polymorphic forms that distinguishes the 
individual See generally National Research Council, H\e Evaluation of Forensic 
5 DNA Evidence (Eds. Pollard et aL, National Academy Press, DC, 1 996). The more 
sites that are analyzed, the lower the probability that the set of polymorphic forms in 
one individual is the same as that in an unrelated individual. Preferably, if multiple 
sites are analyzed, the sites are unlinked. Thus, polymorphisms of the invention are 
often used in conjunction with polymorphisms in distal genes. Preferred 
10 polymorphisms for use in forensics are biallelic because the population frequencies 
of two polymorphic forms can usually be determined with greater accuracy than 
those of multiple polymorphic forms at multi-allelic loci. 

The capacity to identify a distinguishing or unique set of forensic markers in 
an individual is useful for forensic analysis. For example, one can determine 
1 5 whether a blood sample from a suspect matches a blood or other tissue sample from 
a crime scene by determining whether the set of polymorphic forms occupying 
selected polymorphic sites is the same in the suspect and the sample. If the set of 
polymorphic markers does not match between a suspect and a sample, it can be 
concluded (barring experimental error) that the suspect was not the source of the 
20 sample. If the set of markers does match, one can conclude that the DNA from the 
suspect is consistent with that found at the crime scene. If frequencies of the 
polymorphic forms at the loci tested have been determined (e.g., by analysis of a 
suitable population of individuals), one can perform a statistical analysis to 
determine the probability that a match of suspect and crime scene sample would 
25 occur by chance. 

p(ID) is the probability that two random individuals have the same 
polymorphic or allelic form at a given polymorphic site. In biallelic loci, four 
genotypes are possible: AA, AB, BA, and BB. If alleles A and B occur in a haploid 
genome of the organism with frequencies x and y, the probability of each genotype 
30 in a diploid organism is (see WO 95/12607): 

Homozygote: p(AA)= x 2 
Homozygote: p(BB)= y 2 = (1-x) 2 
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Single Heterozygote: p(AB)= p(BA)= xy = x(l-x) 
Both Heterozygotes: p(AB+BA)= 2xy = 2x(l-x) 

The probability of identity at one locus (i.e, the probability that two 
individuals, picked at random from a population will have identical polymorphic 
5 forms at a given locus) is given by the equation: 
p(ID) « (x 2 ) 2 4- (2xy) 2 + (y 2 ) 2 . 

These calculations can be extended for any number of polymorphic forms at 
a given locus. For example, the probability of identity p(ID) for a 3-allele system 
where the alleles have the frequencies in the population of x, y and z, respectively, is 
1 0 equal to the sum of the squares of the genotype frequencies: 
p(ID) = x 4 + (2xy) 2 + (2yz) 2 + (2xz) 2 + z 4 + y 4 

In a locus of n alleles, the appropriate binomial expansion is used to calculate 
p(ID) and p(exc). 

The cumulative probability of identity (cum p(ID)) for each of multiple 
15 unlinked loci is determined by multiplying the probabilities provided by each locus, 
cum p(ID) = p(DDl)p(ID2)p(ID3).... p(IDn) 

The cumulative probability of non-identity for n loci (i.e. the probability that 
two random individuals will be different at 1 or more loci) is given by the equation: 

cum p(nonlD) ~ 1 -cum p(ID). 
20 If several polymorphic loci are tested, the cumulative probability of non- 

identity for random individuals becomes very high (e.g., one billion to one). Such 
probabilities can be taken into account together with other evidence in determining 
the guilt or innocence of the suspect. 

B. Paternity Testing 

25 The object of paternity testing is usually to determine whether a male is the 

father of a child. In most cases, the mother of the child is known and thus, the 
mother's contribution to the child's genotype can be traced. Paternity testing 
investigates whether the part of the child's genotype not attributable to the mother is 
consistent with that of the putative father. Paternity testing can be performed by 

30 analyzing sets of polymorphisms in the putative father and the child. 
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If the set of polymorphisms in the child attributable to the father does not 
match the set of polymorphisms of the putative father, it can be concluded, barring 
experimental error, that the putative father is not the real father. If the set of 
polymorphisms in the child attributable to the father does match the set of 
5 polymorphisms of the putative father, a statistical calculation can be performed to 
determine the probability of coincidental match. 

The probability of parentage exclusion (representing the probability that a 
random male will have a polymorphic form at a given polymorphic site that makes 
him incompatible as the father) is given by the equation (see WO 95/12607): 
1 0 p(exc) = xy(l-xy) 

where x and y are the population frequencies of alleles A and B of a biallelic 
polymorphic site. 

(At a triallelic site p(exc) = xy(l-xy) + yz(l- yz) + xz(l-xz)+ 3xyz(l-xyz))) ) 
where x, y and z and the respective population frequencies of alleles A, B and C). 
1 5 The probability of non-exclusion is 

p(non-exc) = 1 -p(exc) 

The cumulative probability of non-exclusion (representing the value obtained 
when n loci are used) is thus: 

cum p(non-exc) = p(non-excl)p(non-exc2)p(non-exc3).... p(non-excn) 
20 The cumulative probability of exclusion for n loci (representing the 

probability that a random male will be excluded) 
cum p(exc) = 1 - cum p(non-exc). 

If several polymorphic loci are included in the analysis, the cumulative 
probability of exclusion of a random male is very high. This probability can be 
25 taken into account in assessing the liability of a putative father whose polymorphic 
marker set matches the child's polymorphic marker set attributable to his/her father. 

C. Correlation of Polymorphisms with Phenotypic Traits 
The polymorphisms of the invention may contribute to the phenotype of an 
organism in different ways. Some polymorphisms occur within a protein coding 
30 sequence and contribute to phenotype by affecting protein structure. The effect may 
be neutral, beneficial or detrimental, or both beneficial and detrimental, depending 
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on the circumstances. For example, a heterozygous sickle cell mutation confers 
resistance to malaria, but a homozygous sickle cell mutation is usually lethal. Other 
polymorphisms occur in noncoding regions but may exert phenotypic effects 
indirectly via influence on replication, transcription, and translation. A single 
5 polymorphism may affect more than one phenotypic trait. Likewise, a single 
phenotypic trait may be affected by polymorphisms in different genes. Further, 
some polymorphisms predispose an individual to a distinct mutation that is causally 
related to a certain phenotype. 

Phenotypic traits include diseases that have known but hitherto unmapped 
10 genetic components. Phenotypic traits also include symptoms of, or susceptibility 
to, multifactorial diseases of which a component is or may be genetic, such as 
autoimmune diseases, inflammation, cancer, diseases of the nervous system, and 
infection by pathogenic microorganisms. Some examples of diseases which can be 
treated or diagnosed as described herein include, but are not limited to, 

1 5 bradyarrhythmias, tachyarrhythmias, heart failure, such as congestive heart failure, 
congenital heart disease, rheumatic fever, valvular heart disease, cardiomyopathies, 
myocarditides, pericardial diseases, cardiac tumors, cardiac manifestations of 
systemic diseases, and traumatic cardiac injury. Other disorders include 
atherosclerosis, acute myocardial infarction, ischemic heart disease, hypertensive 

20 vascular disease, disorders of the aorta, vascular diseases of the extremities, vessel 
wall disorders, such as various forms of thrombocytopenia, von Willebrand's disease 
and drug-induced platelet dysfunction, and homeostatic disorders relating to vessel 
disease and associated bleeding. Also suitable are thrombotic thrombocytopenic 
purpura, hemolytic-uremic syndrome, Henoch-SchGnlein purpura, capillary fragility, 

25 vascular purpura, metabolic and inflammatory disorders, such as those induced by 
rickettsiae and certain drugs, such as sulfonamides, aortic aneurysm, aortic 
dissection, aortic occlusion, aortitis, atherosclerosis, coronary artery disease, angina, 
myocardial infarction, thrombosis, hemostatic and coagulation disorders, 
hypertension and hypotension. Other disorders include transplant accelerated 

30 vascular restenosis following balloon angioplasty, Raynaud's disease and 
acrocyanosis. 
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Additional disorders include, but are not limited to, disorders of 
neurodegeneration characterized by astrocyte hypertrophy including gliosis, Pick's 
disease, aceroplasminemia, portal-systemic encephalopathy, frontal lobe dementia 
and inherited and acquired ataxias, neurodegenerative diseases of other etiology 
5 including progressive supranuclear palsy, primary progressive aphasia, cortical basal 
degeneration, Alzheimer's disease, Huntington's disease, and Parkinson's disease, 
retinitis pigmentosa and amyotrophic lateral sclerosis. Other disorders include 
epilepsy, stroke, defects of neural migration and differentiation, including Miller- 
Dieker lissencephaly syndrome, and cancer of the brain including astrocytomas and 
10 gliomas, as well as psychological disorders such as schizophrenia. 

Phenotypic traits also include characteristics such as longevity, appearance 
(e.g., baldness, obesity), strength, speed, endurance, fertility, and susceptibility or 
receptivity to particular drugs or therapeutic treatments. 

The correlation of one or more polymorphisms with phenotypic traits can be 
1 5 facilitated by knowledge of the gene product of the wild type (reference) gene. The 
genes in which cSNPs of the present invention have been identified are genes which 
have been previously sequenced and characterized in one of their allelic forms. 

Correlation is performed for a population of individuals who have been 
tested for the presence or absence of a phenotypic trait of interest and for 
20 polymorphic markers sets. To perform such analysis, the presence or absence of a 
set of polymorphisms (i.e. a polymorphic set) is determined for a set of the 
individuals, some of whom exhibit a particular trait, and some of which exhibit lack 
of the trait The alleles of each polymorphism of the set are then reviewed to 
determine whether the presence or absence of a particular allele is associated with 
25 the trait of interest. Correlation can be performed by standard statistical methods 
such as a K-squared test and statistically significant correlations between 
polymorphic form(s) and phenotypic characteristics are noted. For example, it 
might be found that the presence of allele Al at polymorphism A correlates with 
heart disease. As a further example, it might be found that the combined presence of 
30 allele Al at polymorphism A and allele Bl at polymorphism B correlates with 
increased milk production of a farm animal. 
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Such correlations can be exploited in several ways. In the case of a strong 
correlation between a set of one or more polymorphic forms and a disease for which 
treatment is available, detection of the polymorphic form set in a human or animal 
patient may justify immediate administration of treatment, or at least the institution 
5 of regular monitoring of the patient. Detection of a polymorphic form correlated 
with serious disease in a couple contemplating a family may also be valuable to the 
couple in their reproductive decisions. For example, the female partner might elect 
to undergo in vitro fertilization to avoid the possibility of transmitting such a 
polymorphism from her husband to her offspring. In the case of a weaker, but still 

10 statistically significant correlation between a polymorphic set and human disease, 
immediate therapeutic intervention or monitoring may not be justified. 
Nevertheless, the patient can be motivated to begin simple life-style changes (e.g., 
diet, exercise) that can be accomplished at little cost to the patient but confer 
potential benefits in reducing the risk of conditions to which the patient may have 

1 5 increased susceptibility by virtue of variant alleles. Identification of a polymorphic 
set in a patient correlated with enhanced receptiveness to one of several treatment 
regimes for a disease indicates that this treatment regime should be followed. 

For animals and plants, correlations between characteristics and phenotype 
are useful for breeding for desired characteristics. For example, Beitz et aL, US 

20 5,292,639 discuss use of bovine mitochondrial polymorphisms in a breeding 

program to improve milk production in cows. To evaluate the effect of mtDNA D- 
loop sequence polymorphism on milk production, each cow was assigned a value of 
1 if variant or 0 if wildtype with respect to a prototypical mitochondrial DNA 
sequence at each of 17 locations considered. Each production trait was analyzed 

25 individually with the following animal model: 

Y^Ji + YS l + P j + X k + P I + ...p l7 + PE B + a B +c i> 
where Y iJknp is the milk, fat, fat percentage, SNF, SNF percentage, energy 
concentration, or lactation energy record; p. is an overall mean; YS f is the effect 
common to all cows calving in year-season; X k is the effect common to cows in 
30 either the high or average selection line; P, to P 17 are the binomial regressions of 
production record on mtDNA D-loop sequence polymorphisms; PE n is permanent 
environmental effect common to all records of cow n; a„ is effect of animal n and is 
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composed of the additive genetic contribution of sire and dam breeding values and a 
Mendelian sampling effect; and e p is a random residual. It was found that eleven of 
seventeen polymorphisms tested influenced at least one production trait. Bovines 
having the best polymorphic forms for milk production at these eleven loci are used 
5 as parents for breeding the next generation of the herd. 

D. Genetic Mapping of Phenotypic Traits 

The previous section concerns identifying correlations between phenotypic 
traits and polymorphisms that directly or indirectly contribute to those traits. The 
present section describes identification of a physical linkage between a genetic locus 
10 associated with a trait of interest and polymorphic markers that are not associated 
with the trait, but are in physical proximity with the genetic locus responsible for the 
trait and co-segregate with it. Such analysis is useful for mapping a genetic locus 
associated with a phenotypic trait to a chromosomal position, and thereby cloning 
gene(s) responsible for the trait. See Lander et aL % Proc. Natl. Acad. Set (USA) 83, 

15 7353-7357 (1986); Lander et al. y Proc. Natl. Acad. Sci. (USA) 84, 2363-2367 
(1987); Donis-Keller et al, Cell 51, 319-337 (1987); Lander et al. 9 Genetics 121, 
185-199 (1989)). Genes localized by linkage can be cloned by a process known as 
directional cloning. See Wainwright, Med. J. Australia 159, 170-174 (1993); 
Collins, Nature Genetics 1, 3-6 (1992). 

20 Linkage studies are typically performed on members of a family. Available 

members of the family are characterized for the presence or absence of a phenotypic 
trait and for a set of polymorphic markers. The distribution of polymorphic markers 
in an informative meiosis is then analyzed to determine which polymorphic markers 
co-segregate with a phenotypic trait. See, e.g., Kerem et al. 9 Science 245, 1 073-1080 

25 (1989); Monaco et al. 9 Nature 316, 842 (1985); Yamoka et at., Neurology 40, 222- 
226 (1990); Rossiter et al. y FASEB Journal 5, 21-27 (1991). 

Linkage is analyzed by calculation of LOD (log of the odds) values. A lod 
value is the relative likelihood of obtaining observed segregation data for a marker 
and a genetic locus when the two are located at a recombination fraction 0, versus 

30 the situation in which the two are not linked, and thus segregating independently 
(Thompson & Thompson, Genetics in Medicine (5th ed, W.B. Saunders Company, 
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Philadelphia, 1991); Strachan, "Mapping the human genome" in The Human 
Genome (BIOS Scientific Publishers Ltd, Oxford), Chapter 4). A series of 
likelihood ratios are calculated at various recombination fractions (9), ranging from 
0 - 0.0 (coincident loci) to 8 = 0.50 (unlinked). Thus, the likelihood at a given 
5 value of 8 is: probability of data if loci linked at 6 to probability of data if loci 
unlinked. The computed likelihoods are usually expressed as the log 10 of this ratio 
(i.e., a lod score). For example, a lod score of 3 indicates 1000:1 odds against an 
apparent observed linkage being a coincidence. The use of logarithms allows data 
collected from different families to be combined by simple addition. Computer 

10 programs are available for the calculation of lod scores for differing values of 8 

(e.g., LIPED, MLINK (Lathrop, Proc. Nat. Acad. Sci. (USA) 81, 3443-3446 (1984)). 
For any particular lod score, a recombination fraction may be determined from 
mathematical tables. See Smith et aL, Mathematical tables for research workers in 
human genetics (Churchill, London, 1961); Smith, Ann. Hum. Genet. 32, 127-150 

15 (1968). The value of 0 at which the lod score is the highest is considered to be the 
best estimate of the recombination fraction. 

Positive lod score values suggest that the two loci are linked, whereas 
negative values suggest that linkage is less likely (at that value of 0) than the 
possibility that the two loci are unlinked. By convention, a combined lod score of 

20 +3 or greater (equivalent to greater than 1000: 1 odds in favor of linkage) is 

considered definitive evidence that two loci are linked. Similarly, by convention, a 
negative lod score of -2 or less is taken as definitive evidence against linkage of the 
two loci being compared. Negative linkage data are useful in excluding a 
chromosome or a segment thereof from consideration. The search focuses on the 

25 remaining non-excluded chromosomal locations. 

III. Modified Polypeptides and Gene Sequences 

The invention further provides variant forms of nucleic acids and 
corresponding proteins. The nucleic acids comprise one of the sequences described 
in Figures 5A-5QQQQQQQ, column 1 1, in which the polymorphic position is 
30 occupied by one of the alternative bases for that position. Some nucleic acids 
encode full-length variant forms of proteins. Similarly, variant proteins have the 
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prototypical amino acid sequences encoded by nucleic acid sequences shown in 
Figures 5A-5QQQQQQQ, column 11, (read so as to be in-frame with the full-length 
coding sequence of which it is a component) except at an amino acid encoded by a 
codon including one of the polymorphic positions shown in Figures 5A- 
5 5QQQQQQQ. That position is occupied by the amino acid coded by the 
corresponding codon in any of the alternative forms shown in Figures 5A- 
5QQQQQQQ. 

Variant genes can be expressed in an expression vector in which a variant 
gene is operably linked to a native or other promoter. Usually, the promoter is a 

10 eukaryotic promoter for expression in a mammalian cell. The transcription 

regulation sequences typically include a heterologous promoter and optionally an 
enhancer which is recognized by the host. The selection of an appropriate promoter, 
for example tip, lac, phage promoters, glycolytic enzyme promoters and tRNA 
promoters, depends on the host selected. Commercially available expression vectors 

15 can be used. Vectors can include host-recognized replication systems, amplifiable 
genes, selectable markers, host sequences useful for insertion into the host genome, 
and the like. 

The means of introducing the expression construct into a host cell varies 
depending upon the particular construction and the target host. Suitable means 

20 include fusion, conjugation, transfection, transduction, electroporation or injection, 
as described in Sambrook, supra. A wide variety of host cells can be employed for 
expression of the variant gene, both prokaryotic and eukaryotic. Suitable host cells 
include bacteria such as E. coli, yeast, filamentous fungi, insect cells, mammalian 
cells, typically immortalized, e.g., mouse, CHO, human and monkey cell lines and 

25 derivatives thereof. Preferred host cells are able to process the variant gene product 
to produce an appropriate mature polypeptide. Processing includes glycosylation, 
ubiquitination, disulfide bond formation, general post-translational modification, and 
the like. 

The protein may be isolated by conventional means of protein biochemistry 
30 and purification to obtain a substantially pure product, i.e., 80, 95 or 99% free of cell 
component contaminants, as described in Jacoby, Methods in Enzymology Volume 
104, Academic Press, New York (1984); Scopes, Protein Purification, Principles 



WO 00/58519 



-32- 



PCT/US00/08440 



and Practice, 2nd Edition, Springer- Verlag, New York (1987); and Deutscher (ed), 
Guide to Protein Purification, Methods in Enzymology, Vol, 182 (1990). If the 
protein is secreted, it can be isolated from the supernatant in which the host cell is 
grown. If not secreted, the protein can be isolated from a lysate of the host cells. 
5 The invention further provides transgenic nonhuman animals capable of 

expressing an exogenous variant gene and/or having one or both alleles of an 
endogenous variant gene inactivated. Expression of an exogenous variant gene is 
usually achieved by operably linking the gene to a promoter and optionally an 
enhancer, and microinjecting the construct into a zygote. See Hogan et aL t 

1 0 "Manipulating the Mouse Embryo, A Laboratory Manual," Cold Spring Harbor 

Laboratory. Inactivation of endogenous variant genes can be achieved by forming a 
transgene in which a cloned variant gene is inactivated by insertion of a positive 
selection marker. See Capecchi, Science 244, 1288-1292 (1989). The transgene is 
then introduced into an embryonic stem cell, where it undergoes homologous 

1 5 recombination with an endogenous variant gene. Mice and other rodents are 
preferred animals. Such animals provide useful drug screening systems. 

In addition to substantially full-length polypeptides expressed by variant 
genes, the present invention includes biologically active fragments of the 
polypeptides, or analogs thereof, including organic molecules which simulate the 

20 interactions of the peptides. Biologically active fragments include any portion of the 
full-length polypeptide which confers a biological function on the variant gene 
product, including ligand binding, and antibody binding. Ligand binding includes 
binding by nucleic acids, proteins or polypeptides, small biologically active 
molecules, or large cellular structures. 

25 Polyclonal and/or monoclonal antibodies that specifically bind to variant gene 

products but not to corresponding prototypical gene products are also provided. 
Antibodies can be made by injecting mice or other animals with the variant gene 
product or synthetic peptide fragments thereof. Monoclonal antibodies are screened 
as are described, for example, in Harlow & Lane, Antibodies, A Laboratory Manual, 

30 Cold Spring Harbor Press, New York (1988); Goding, Monoclonal antibodies, 
Principles and Practice (2d ed.) Academic Press, New York (1986). Monoclonal 
antibodies are tested for specific immunoreactivity with a variant gene product and 
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lack of immunoreactivity to the corresponding prototypical gene product. These 
antibodies are useful in diagnostic assays for detection of the variant form, or as an 
active ingredient in a pharmaceutical composition. 

IV. Kits 

5 The invention further provides kits comprising at least one allele-specific 

oligonucleotide as described above. Often, the kits contain one or more pairs of 
allele-specific oligonucleotides hybridizing to different forms of a polymorphism. 
In some kits, the allele-specific oligonucleotides are provided immobilized to a 
substrate. For example, the same substrate can comprise allele-specific 

10 oligonucleotide probes for detecting at least 10, 100 or all of the polymorphisms 
shown in Figures 5A-5QQQQQQQ. Optional additional components of the kit 
include, for example, restriction enzymes, reverse-transcriptase or polymerase, the 
substrate nucleoside triphosphates, means used to label (for example, an avidin- 
enzyme conjugate and enzyme substrate and chromogen if the label is biotin), and 

1 5 the appropriate buffers for reverse transcription, PCR, or hybridization reactions. 
Usually, the kit also contains instructions for carrying out the methods. 

From the foregoing, it is apparent that the invention includes a number of 
general uses that can be expressed concisely as follows. The invention provides for 
the use of any of the nucleic acid segments described above in the diagnosis or 

20 monitoring of diseases, such as coronary artery disease, diabetes, coagulation 

disorders, lipid metabolism disorders, energy metabolism disorders, diseases of the 
blood, blood vessels and cardiovascular system, and infection by microorganisms, as 
well as psychological disorders (e.g., bipolar disorder, psychiophrenia). The 
invention further provides for the use of any of the nucleic acid segments in the 

25 manufacture of a medicament for the treatment or prophylaxis of such diseases. The 
invention further provides for the use of any of the DNA segments as a 
pharmaceutical. 



30 



While this invention has been particularly shown and described with 
references to preferred embodiments thereof, it will be understood by those skilled 
in the art that various changes in form and details may be made therein without 
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CLAIMS 

What is claimed is: 

1 . A nucleic acid molecule selected from the group consisting of the genes listed 
in Figures 5A-5QQQQQQQ, wherein said nucleic acid molecule is at least 5 

5 nucleotides in length and comprises a polymorphic site identified in Figures 

5A-5QQQQQQQ, wherein a nucleotide at the polymorphic site is different 
from a nucleotide at the polymorphic site in a corresponding reference allele. 

2. A nucleic acid molecule according to Claim 1, wherein said nucleic acid 
molecule is at least 10 nucleotides in length. 

10 3. A nucleic acid molecule according to Claim 1, wherein said nucleic acid 
molecule is at least 20 nucleotides in length. 

4. A nucleic acid molecule according to Claim 1 , wherein the nucleotide at the 
polymorphic site is the variant nucleotide for the gene listed in Figures 5A- 
5QQQQQQQ. 

15 5. An allele-specific oligonucleotide that hybridizes to a portion of a gene 
selected from the group consisting of the genes listed in Figures 5A- 
5QQQQQQQ, wherein said portion is at least 5 nucleotides in length and 
comprises a polymorphic site identified in Figures 5A-5QQQQQQQ, wherein 
a nucleotide at the polymorphic site is different from a nucleotide at the 

10 polymorphic site in a corresponding reference allele. 

6. An allele-specific oligonucleotide according to Claim 5 that is a probe. 

7. An allele-specific oligonucleotide according to Claim 5, wherein a central 
position of the probe aligns with the polymorphic site of the portion. 



WO 00/58519 



-36- 



PCT/USOO/08440 



8. An allele-specific oligonucleotide according to Claim 5 that is a primer. 

9. An allele-specific oligonucleotide according to Claim 8, wherein the 3* end of 
the primer aligns with the polymorphic site of the portion. 

10. An isolated gene product encoded by a nucleic acid molecule according to 
5 Claim 1. 

11. A method of analyzing a nucleic acid sample, comprising obtaining the 
nucleic acid from an individual sample; and determining a base occupying any 
one of the polymorphic sites shown in Figures 5A-5QQQQQQQ. 

12. A method according to Claim 1 1, wherein the nucleic acid sample is obtained 
1 0 from a plurality of individuals, and a base occupying one of the polymorphic 

positions is determined in each of the individuals, and the method further 
comprising testing each individual for the presence of a disease phenotype, 
and correlating the presence of the disease phenotype with the base. 
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1001 


1002 


Assay Sequence 


CACATCCATACTGCCTGAGTcagccccgggttacgcctgttgtcccggcataaccattgctagca 

caccctttccctctcagaagtgccccggtttgaatgaaacctcttcgtgatccccttggaggtca 

actctgagggacccagaaactgccttttgactgcatttagtactccatgaagtcaccctcatttc 

tttttcattccaggtgcatagcgtaatgtccatgttgttctacactctgatcacagcttttctga 

tcggcatacaggc(g/algaaccacactcagagagcaatgtccctgcaggacacaccatccccca 

agtccactggactaaacttcagcattcccttgacactgcccttcgcagagcccgcagcgccccgg 

cagcggcgatagctgcacgcgtggcggggcagacccgcaacattactgtggaccccaggctgttt 

aaaaagcggcgactccgttcaccccgtgtgctgtttagcacccagcctccccgtgaagctgcaga 

cactcaggatctggacttcgaggtcggtggtgctgcccccttcaacaggactcacaggagcaagc 

ggtcatcatcccatcccatcttccacaggggcgaattctcggtgtgtgacagtgtcagcgtgtgg 

gttggggataagaccaccgccacagacatcaagggcaaggaggtgatggtgctgggagaggtgaa 

cattaacaacagtgtattcaaacagtacttttttgagaccaagtgccgggacccaaatcccgttg 

acagcgggtgccggggcattgactcaaagcactggaactcatattgtaccacgactcacaccttt 

gtcaaggcgctgaccatggatggcaagcaggctgcctggcggtttatccggatagatacggcctg 

tgtgtgtgtgctcagcaggaaggctgtgagaagagcctgacctgccgacacgctccctcccccta 
cccc t tc taCACTCTCCTGGGCCCCe 


mACAGAGGAGCTGACGTTTGctacacatctacaagtatgcataggagctccgcggaggccagt 

gagaggccctccaggagcagaactaattccacaattacttgaccaagttgggggattatttgtgg 

glg/cl taactgcagtgcagtatggagtcctcttggggacagttagagccataccatttgatcta 

tagtcacataagaacaaacaataaaaagaaaagacatgcttaagagtgaaagagaaagggaggga 

gaaaagaaggaagggtggatggaaggacactagcttagtaaggggtcaactttggattctatttc 

tggttcagttttcatttgtgacttcagtgctttagtgttagttcattttctgtgaatcagtttcc 

ttatatgtgaaataaatatgataaatcctaattgaactcacagactcatgagaagatagaagtga 

acacattttaaaaacatcacacaaagaggaactattatgtggtccacatttatatatgtggggta 

gcgtctgaagaggtgcctggactaagatggtcccagagccacaaggtttttgccaaacatgacgc 

tttgtgaattcataacaagggctccaagtcaccagatcttagagctgacccagtgcactgtctga 

aagggggtaccagttctgaggcttcaagacatgtccccagcagatcttccccgtgccttcccaga 

ggattcaaaactgttgagcaggacggcaccatcacatcaaggcacaagtgccagggagaggtgtt 

aaactctccccaccaacctccctggtacacacatggacacttaccacctccctcagccgccttaa 

gcttcagagaactcaaaggactctgtaagtgatgtctccaagctcatatcgaactactgggcaaa 

atttcaggggctctgtcacttcctggagaagctcggatggggtgaccacacatccatactgcctg 
ag t cagcccCGGGTT ACGCCTGTTGTC 


Reverse 
Primer 
iS* -> 3*) 


AGGGGCCCAGG 
AGAGTG 
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GACAACAGGCG 
TAACCCG 
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Forward 
Primer 
(5' -> 3* 


CACATCCATJ 
TGCCTGAGT 




TTTACAGAGG 
GCTGACGTTT 
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U jJ « 
U y CO 
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y JJ CQ 
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O iJ D 
<0 i-» u 
O U JJ 
U ID U 

o y o 
y y u 
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D4J CD 

n to io 
U (S u 
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Jj Jj Jj 
Oiu u 
J D Ol 
CO Ol (Q 

<o o> o 

U JJ JJ 
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JJ JJ JJ 
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01 U <D 

^ co y 
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Ol O Ol 
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U (0 O) 
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u cn jj < 
y oi u p 
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jJ jJ U t-* 

o o» « 

« JJ u CJ 
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Ol U JJ 
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